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Abstract. Measurements of C0 2 , N 2 0 and CH 4 are analyzed to 
define hemispheric average vertical exchange rates in the lower 
stratosphere from November 1992 to October 1993. Effective 
vertical diffusion coefficients were small in summer, < 1 m 2 s _1 at 
altitudes below 25 km; values were similar near the tropopause in 
winter, but increased markedly with altitude. The analysis 
suggests possibly longer residence times for exhaust from 
stratospheric aircraft, and more efficient transport from 20 km to 
the middle stratosphere, than predicted by many current models. 
Seasonally-resolved measurements of stratospheric C0 2 and N 2 0 
provide significant new constraints on rates for global-scale 
vertical transport. 


one tracer concentration against another give compact curvilinear 
relationships, often uniform on global scales [e.g., Fahey et al., 
1990]. Plumb and Ko [1992] argued that compact global relation- 
ships arise by advection of trace species in the global residual cir- 
culation, modified by quasi-isentropic mixing by planetary -scale 
waves that propagate from the troposphere in winter [cf. Holton, 
1986, Mahlman et al., 1986; Plumb and Mahlman, 1987; Schoe- 
berl and Hartmann, 1991]. If quasi-isentropic mixing were much 
faster than non-conservative processes, isopleths for long-lived 
tracers would be parallel and transport of one tracer relative to 
another could be represented as a flux (F z *) in one dimension, nor- 
mal to the isopleths, 

F z . =- K(Z*) , where (la) 


Introduction 

Future high speed civil transports (HSCTs) are expected to 
inject H 2 0 and nitrogen oxides into the stratosphere. The 
residence time for the exhaust in the stratosphere will determine 
the magnitude of the associated perturbation to atmospheric chem- 
istry. Residence times for stratospheric gases are currently 
estimated using observations of radioactive and volcanic debris, 
mostly in particulate form. These data reflect stratospheric tran- 
sport rates only for particular single-point injections in the past. 

To predict the impact of HSCTs, we must understand rates for 
meridional dispersion and vertical transport of exhaust continu- 
ously emitted over the globe. Boering et al. [fAij issue ] docu- 
mented the propagation through the stratosphere of seasonal 
cycles and long-term trends in tropospheric C0 2 [Keeling et al., 
1989; Conway et al.,1988], relative to N 2 0. Here we use data for 
N 2 0, C0 2 , and CH 4 obtained from the ER-2 during SPADE to 
derive seasonally-resolved rates for vertical exchange in 1992-93, 
averaged over midlatitudes in the Northern hemisphere strato- 
sphere, and we discuss implications for predicted impacts of 
HSCTs. 

Conceptual framework 

Concentrations of stratospheric trace gases fluctuate with lati- 
tude, altitude, and time. However, variations for long-lived 
species are usually correlated [Ehhalt et al., 1983]: scatterplots of 
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Here a is the tracer mole fraction, % the stream function, and K a 
mixing tensor along (y,z), the meridional and vertical directions, 
respectively. Brackets denote the global average on an isopleth 
and Z* (= 71n( 1000/p) (km), p=isopleth pressure (mb)) is the iso- 
pleth pressure altitude at a specified reference latitude . 

The data conform closely to this idealized model. Scatterplots 
of C0 2 [Boering et al. this wsk*] and CH 4 vs. N 2 0 were remark- 
ably compact. Seasonal and interannual variations of C0 2 pro- 
pagated uniformly upward from the troposphere (seasonal changes 
for CH 4 were undetectable when compared with earlier measure- 
ments [see Schauffler et al., 1993]). Hence we use Eq. (1) with 
data for C0 2 , CH^ and N 2 0 to derive mean rates for vertical tran- 
sport in the stratosphere over the domain of the observations (15- 
60°N), analogous to early "eddy-diffusion" models [Lettau, 1951; 
Colegrove et al., 1965; Wofsy and McElroy, 1973; Hunten, 1975]. 
Here the C0 2 /CH4-N 2 0 coordinate system eliminates effects of 
variance due to reversible displacements of the tracer fields 
[Ehhalt et al., 1983], allowing us to exploit these displacements to 
infer tracer distributions above flight altitudes. 


Mean K(Z*) for 1992-93 


We derive K(Z*) by integrating the continuity equation. 
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Here C = CH 4 + C0 2 , and the last term (with H=6 km) approxi- 
mates the flux at Z*=Z* top . Eq. 3 was evaluated using data for 
AC0 2 /AN 2 0 and AC0 2 (N 2 0)/At [Fig. 2, Boering et al.Jhis issue], 
and ACH 4 /AN 2 0 = 0.0043 ppm/ppb. A profile for N 2 0(Z*) was 
adopted from data for 24±4° (Fig. lb), the reference latitude 
chosen to approximate global mean photolysis rates for N 2 0. The 
profile for N 2 0(Z*) from Fig. lb was used to derive K(Z*) for 42 
> Z* > 30 km, the interval with 70% of the N 2 0 gradient, using 
Eq. 3 with dC(N 2 0(Z*))/dt replaced by A from Eq. (4) below. 
Including CH 4 increases vertical gradients by -15%, with little 
effect on time derivatives, hence, neglect of C0 2 from CH 4 oxida- 
tion would give K(Z*) about 15% too large. 

Profiles of K(Z*) (Fig. la) show a bottleneck just above the 
tropopause, implied by slow propagation of the seasonal oscilla- 
tion. The fall minima and spring maxima for C0 2 require 4-7 
months to propagate 2-3 km, i.e. from the tropopause to N 2 0=300 
ppb. The bottleneck was inferred also in early studies of transport 
of heat, nuclear bomb debris, CH 4 and N 2 0 [e.g. Holton, 1986; 
Wofsy and McElroy, 1973; Hun ten, 1975], from steep gradients 
observed just above the tropopause. The increase of K(Z*) above 
the tropopause in winter is implied by the rise in C0 2 for N 2 0 < 
200 ppb; very slow vertical exchange at these levels is indicated in 
summer by virtually unchanged C0 2 from May to November. 

The composite profile adopted for N 2 0 is not constrained to 
be consistent with the C0 2 -N 2 0 relationships in SPADE. Figure 
2a shows consistent solutions for C0 2 (Z*) from Eq. (2) and for 
N 2 0 from 



Figure la. Vertical profiles of K(Z*): summer (open □) and 
winter (filled □) using Eq. 3 and data for C0 2 , CH 4 , and N 2 0 
(Z*<26 km); solid line, using satellite data for N 2 0(Z*) (Z* > 30 
km); dotted line, interpolated values. 

Figure lb. Vertical profile of N 2 0 for 15-25N. The line is a fit 
(locally-weighted least squares, "lowess") to data: Z* >30 km, 
satellite observations November, 1992 [CLAES/UARS, A. Roche 
and J. Mergenthaler, private communication, 1994; 
ATMOS/ATLAS, M. Gunson, private communication, 1994]; 30 
> Z* > 20 km, tropical balloon data [Goldan et al., 1981]; Z* < 20 
km, ER-2, 20-28N, on dates indicated. Dotted lines show N 2 0 
computed from the 1-D model, using K(Z*) from Eq. (3) and with 
0.5x this valueC'O.SxK^”). 
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Figure 2a. Relationships for C0 2 vs N 2 0 computed using the 1- 
D model, with K(Z*) from Fig. la (sK^), and scaled by 0.5, for 
November 1992 (upper) and May 1993 (lower). 

Figure 2b. Concentrations of C0 2 for given N z O from the 1 -D 
model with K(Z*)=0.55*K dir (lines), compared to observations 
(symbols). Data for C0 2 at the lowest altitude unambiguously 
above the tropopause (□ "L-S") [Boering et al, this issue] and the 
boundary condition for the 1-D model ( — , M Surf', mean of 


NOAA/CMDL data at Samoa and Mauna Loa [T. Conway, private 
communication, 1994] delayed by 1 month) are also shown. 
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Here we set tropospheric N 2 0 to 325 ppb, as given by the ER-2 
ATLAS instrument. We set C0 2 at the bottom boundary to the 
mean of surface data from Mauna Loa and Samoa [T. Conway, 
private communication, 1994], delayed by one month to allow for 
tropospheric mixing, which approximates C0 2 observations just 
above the tropopause (Fig. 2b). In Eq. (4), A = J N20 + k 5 [0( J D)] 
for 15N at equinox, J N2 q is the photolysis rate (sec -1 ), and k 5 is 
the rate constant for 
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0( 1 D) + N 2 0->2N0 (5a) 

0( l D) + N 2 0^>N 2 + 0 2 . (5b) 


The calculation commenced in 1982. The initial C0 2 -N 2 0 rela- 
tionship was specified by subtracting 15 ppm C0 2 (the global 
mean increase in the decade) from the observations for November 
1992. Winter and summer profiles for K(Z*) dk (Fig. la), or 
0.5xK(Z*) dir , were repeated each year. 

Observed C0 2 -N 2 0 relationships are bracketed by models 
using lx and 0.5xK(Z*) dir (Fig. 2a). The 1-D model simulates 
observed propagation of the C0 2 signal through the lower strato- 
sphere remarkably well (using N 2 0 as the vertical coordinate), 
with the best fit for K(Z*)=0.55xK(Z*) dif (Fig. 2b), and it projects 
changes observed in early 1994. Details of N 2 0-C0 2 variations 
near the tropopause (290<N 2 O<325 ppb) are not clearly del- 
ineated: these data were obtained on ascent and descent, adversely 
affecting the diode laser spectrometer. 

Values of N 2 0 are somewhat lower than observed for the 
profile of K(Z*) giving the best fit to C0 2 -N 2 0 data (Fig. lb). 
This difference likely reflects errors in the specified vertical profile 
for N 2 0, leading to an inaccurate global mean vertical flux of 
N 2 0. The discrepancy contributes uncertainty to estimates of the 
rate of upward spread of HSCT exhaust; it could be resolved by 
extending measurements of C0 2 , CH 4 , and N 2 0 to the middle 
stratosphere. 

A recent 3-D model [Hall and Prather, 1993] simulated the 
mean difference between stratospheric and tropospheric C0 2 , 
however, C0 2 evolved uniformly vs N 2 0 through the year, lack- 
ing the observed seasonal asymmetry. The model gave qualita- 
tively correct gradients near the tropopause, but finer vertical reso- 
lution (=0.5 km, vs 3 km in the model) is needed for detailed com- 
parison. 

NO y from stratospheric aviation 


The concentration of NO y was computed by solving 

- -3|-pK(Z*)^^=P(N 2 0(Z*)>-p(N0 y ) 2 (6) 

dt p d/- oLr 


where P(N 2 0(Z*)) denotes production of odd nitrogen by (5a) or 
by aircraft, and p(NO y ) 2 is the sink for NO y , 


NO + hv — > N+O 
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averaged along an N 2 0 isopleth. Solutions to (6) reproduce 
observed correlations for N 2 0-NO y (Figure 3a) [e.g., Loewenstein 
et al., 1993], virtually independent of K(Z*). 

Figure 3b shows computed springtime perturbations to NO y 
produced by hypothetical HSCTs injecting NO y at Z*=20 km for 
10 years at a rate of 1.45xl0 8 cm 2 s~\ 60% of total fleet emis- 
sions for the Mach 2.4/EI 15 fleet scenario [Baughcum et al., 
1993] (to approximate northern hemisphere cruise-level inputs). 
Each summer, excess NO y accumulated near 20 km, then spread 
to higher altitudes, or was removed to the troposphere, in winter. 
The maximum ANO y for the best-fit profile (0.5xK(Z*) dif ) was =7 
ppb (Fig. 3b), 60% larger than predicted by most 2-D models. 
Values for ANO y are insensitive to K(Z*) above, and proportional 
to K(Z*) below, the injection altitude (see curve for lxK(Z*) dir , 
Fig. 3b). The perturbation is reduced by 50% for injection just 2 
km lower (al Z*= 18 km). 
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Figure 3a. N 2 0-NO y relationship computed from Eq. (6) and a 
linear regression for N 2 O>120 ppb, NO y = 25.3 - 0.076*N 2 0. 
Figure 3b. ANO y for HSCTs from the 1-D model for KfZ*)^ 
(from Fig. la) and 0.5xK(Z*) dir , assuming injection at Z* = 20 km 
of 60% of the fleet emissions for the nominal Mach 2.4/EI 15 
scenario of Baughcum et al. [1993] (approximating northern hemi- 
sphere input at cruise). Results are shown from the AER 2-D 
model [Ko and Douglass, 1993] for the full scenario, averaged 
over the Northern hemisphere. 

At altitudes higher than injection, ANO y in 2-D models is 
smaller than in 1-D models, because exhaust does not spread uni- 
formly along N 2 0 isopleths, giving rise to a local maximum in 
flight corridors. It is presently not clear if the 1-D model 
represents a global mean, as envision by Plumb and Ko [1992], or 
if the tropics represent a distinct regime essentially uncoupled 
from midlatitudes. Data for C0 2 , CH 4 , N 2 0 and NO y are 
presently scarce or nonexistent for altitudes above 20 km, prevent- 
ing confident prediction of ANO y above flight levels; this is highly 
significant due to the enhanced reactivity of NO x at high altitudes. 

Conclusions 

Temporal variations of tropospheric C0 2 , relative to N 2 0 and 
CH 4 , provide a unique signal for determining rates for transport in 
the lower and middle stratosphere. Seasonally-resolved data for 
C0 2 , N 2 0 and NO y are needed above 20 km to define rates for 
vertical advection and quasi-isentropic mixing, and fine-scale data 
near the tropopause are needed to define rates for transport in the 
key ’’bottleneck” region. Interannual variations, rates for 
interhemispheric exchange, and possible effects of the quasi- 
biennial oscillation need to be investigated. The present analysis 
shows that systematic tracer measurements over a period of years 
promise to help resolve these issues, providing essential con- 
straints on models used to predict impacts of future HSCTs. 
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